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ABSTRACT 

The HST/ACS colour-magnitude diagrams (CMD) of the populous LMC star 
cluster NGC 1751 present both a broad main sequence turn-off and a dual clump of red 
giants. We show that the latter feature is real and associate it to the first appearance 
of electron-degeneracy in the H-exhausted cores of the cluster stars. We then apply to 
the NGC 1751 data the classical method of star formation history (SFH) recovery via 
CMD reconstruction, for different radii corresponding to the cluster centre, the cluster 
outskirts, and the underlying LMC field. The mean SFH derived from the LMC field 
is taken into account during the stage of SFH-recovery in the cluster regions, in a 
novel approach which is shown to significantly improve the quality of the SFH results. 
For the cluster centre, we find a best-fitting solution corresponding to prolonged star 
formation for a for a timespan of 460 Myr, instead of the two peaks separated by 200 
Myr favoured by a previous work based on isochrone fitting. Remarkably, our global 
best-fitting solution provides an excellent fit to the data - with ^ind residuals close 
to the theoretical minimum - reproducing all the CMD features including the dual red 
clump. The results for a larger ring region around the centre indicate even longer star 
formation, but in this case the results are of lower quality, probably because of the 
differential extinction detected in the area. Therefore, the presence of age gradients 
in NGC 1751 could not be probed. Together with our previous findings for the SMC 
cluster NGC 419, the present results for the NGC 1751 centre argue in favour of 
multiple star formation episodes (or continued star formation) being at the origin of 
the multiple main sequence turn-offs in Magellanic Cloud clusters with ages around 
1.5 Gyr. 

Key words: Stars: evolution - Hertzsprung-Russell (HR) and C-M diagrams 



1 INTRODUCTION 

A few stars clusters in the Magellanic Clouds present in 
their CMDs, i n addition to multiple main sequence turn- 
offs (MMSTOilMackev k Brobv Nielsenll2Q07l:lMackev et al 



2008 1: iMilone et al.ll2009l : iGoudfrooii et a l"2009':'Glat t et al 
2008 ), also dual red clumps (jPiatti et al. 1999; Girar di et al 
20091 ). The MMSTO features can be interpreted either as 



the presence of different generations of stars spanning sev- 



Based on observations with the NASA/ESA Hubble Space Tele- 
scope, obtained at the Space Telescope Science Institute, which 
is operated by the Association of Universities for Research in As- 
tronomy, Inc., under NASA contract NAS5-26555 



eral 10^ yr in t hese c l usters ([Mackev Brobv Nielsenll2007l; 
Mackev et all l2008l : iMilone et al.1 l2009l : IGoudfrooii et al.1 
20091 ). or as the manifestation of some other effect in- 
trinsic to coeval stars such as a disper sion in rota- 
tional velocities (iBastian de MinkI l2009l , see however 
iGirardi, Eggenberger, Miglioll201lh . However, the simul- 
taneous presence of a dual red clump feature favours the for - 
mer interpretation ([Girardi et al ] |2009l : iRubele et al1l2010l ) : 
it indicates a modest spread in the core mass of stars leav- 
ing the main sequence which is well compatible with the 
age spread of a few 10^ yr deduced from the shape of the 
MMSTOs. 

The SMC star cluster NGC 419 is presently the most 
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striking example of a cluster containing a dual red clump. 
iRubele et alj (|2010l l demonstrated that the assumption of 
an extended star formation history (SFH) in NGC 419, ex- 
plored by means of the classical method of SFH-recovery via 
CMD-reconstruction, produces indeed a remarkably good 
quantitative description of the observed CMDs. The SFH 
was found to extend over a period of 700 Myr. The same 
analysis has produced quite stringent limits to the cluster's 
distance, reddening, and metallicity. 

Although not explicitly discussed by I Rubele et alJ 

(|2010l ). the dual red clump of NGC 419 has played an impor- 
tant role in limiting the family of stellar models that could be 
fit in the process of CMD-reconstruction, because dual red 
clumps can only happen within a relatively narrow interval 
of ages. In this paper, we examine the case of the LMC sta r 
cluster NGC 1751, which as noted bv lCirardi et al.1 (l2009l ). 
does also appear to present a dual red clump, and hence 
should be a good target for the CMD-reconstruction tech- 
nique. We will use the extremely accurate data available 
from HST /ACS, and analyses tec hniques simila r to th ose 
applied bv lCirardi et al.1 (l2009l ) and lRubele et al.1 (l2010h for 
deriving the SFH. Sect [2] will briefly present the data and 
discuss the reality of the dual red clump. The next sections 
will present the SFH-recovery method and its application to 
the NGC 1751 surrounding LMC field (Sect. [3]) and cluster 
area (Sect.|4]). Sect. [5] draws the final conclusions. 
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Figure 1. Map of the stars used in this work, in the xy plane of 
the ACS/WFC images. The scale is of about 0.05 arcsec/pix. The 
observed stars have been grouped in areas corresponding to the 
LMC field (red) and, for NGC 1751, an inner "Centre" (green) 
and outer "Ring" (blue). 



2 THE NGC 1751 DATA AND ITS DUAL RED 
CLUMP 

2.1 Data and photometry 

The dataset used in this paper comes from GO- 10595 (PI: 
Goudfrooij), and consists of one short and two long expo- 
sures in F435W, F814W, and F555W with smah dither pat- 
tern to avoid the ACS/WFC gap between two WFC chips. 
A detailed descriptio n of the observ ations and photometry 
is given in Goudfr ooii et al.l (|2009h . Nevertheless, in this 
paper we use the simultaneous ePSF fitting technique as 
it described in Anderson et.al. (2008), which fits the PSF 
simultaneously on al l exposures /obser v ation s of the clus- 
ter. Differently from ICoudfrooii et al.l (|2009l l. the Charge 
Transfer Efficiency (CTE) correction was performed using 
iRiess Mac3 (l2004l ) formula (ACS-ISR 2005). The derived 
photomet ry was calibrated i nto th e Vegamag system as de- 
scribed in ICoudfrooii et al.1 (|2009l ). 

Fig. [1] shows a spatial representation of the stars we 
analize in this work. We have initially located the center of 
NGC 1751 at the position Xc = 2200 pix, yc = 3100 pix. 
Based on the appearance of the CMDs at varying radii r 
from this center, we have defined regions of the ACS/WFC 
images corresponding to 

• the LMC Field, for r > 1500 pix (73.5 arcsec); 

• the main body of the NGC 1751 cluster, for r < 800 pix 
(39.2 arcsec), which is further divided into: 

- the Centre, for r < 350 pix (17.2 arcsec); 

- the Ring, for 350 < r < 800 pix (17.2 <r < 39.2 arc- 
sec). 

These regions are depicted in Fig. [T] The Centre, Ring 
and Field regions have areas of 0.385, 1.63, and 10.30 pix^ 
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Figure 2. The logarithm of stellar density as a function of radius 
from the NGC 1751 centre. Error bars are the random errors. 

(5.3, 22.6, 143 arcmin^), respectively. Figure [2] shows how 
the stellar density varies as a function of radius from the 
NGC 1751 centre, taking into account only the stars of 
FSIAW < 22, for which the photometry should be close 
to complete. The figure clearly shows the flattening of the 
density for r > 1500 pix, which indicates that indeed that 
is a good choice for defining the LMC Field. 
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Figure 3. The CMDs for NGC 1751 as derived from the 
ACS/WFC data, using both F435W - F8UW (top panels) and 
F435W - FSUW (bottom panels) colours versus the F814H^ 
magnitude. Panels from left to right present data for the clus- 
ter Centre and Ring, and LMC Field. As a reference to the eye, 
the panels also show the position of 1.12 and 1.42-Gyr isochrones 
of metal content Z = 0.008, shifted by (m-M)o = 18.50 and 
Ay = 0.7, together with the expected location of equal-mass bina- 
ries along the main sequence (continuous lines). The tiny crosses 
at the leftmost extreme of the left panel are la error bars derived 
from artificial star tests in the cluster Centre (see Sect. [221) • The 
error bars for the Ring and Field are not shown in the figure; they 
are of about the same size for the brightest magnitudes, becoming 
just ~ 25 % smaller for the faintest magnitudes. 

In this work, we will analyse both the Centre and Ring 
regions, which present a good density contrast with respect 
to the field. Although the cluster clearly extends up to a 
radius of 1500 pix, the region with 800 < r < 1500 pix will 
not be considered further. 

Fig. [3] shows the ACS data for the different regions of 
NGC 1751, in the F814W vs. F435W - F814W and F814W 
vs. F555W - F814W CMDs. These plots will be used as a 
reference in our analysis. 

The CMDs for the cluster Centre show very clearly the 



broad main sequence turn-off of this cluster, its dual red 
clump, and other well-known CMD features such as the se- 
quence of binaries parallel to the main sequence, and the 
RGB, subgiants, and early-AGB bump. A simple compar- 
ison between the CMDs for the Centre and Field reveals 
that the field contamination in the cluster central regions 
is close to negligible. This is clear already looking at the 
star counts in the red clump: the 7.15 arcmin^ of Field 
contain 189 red clump stars (here defined as stars with 
18.05 < FSUW < 19.15, F435W - F8UW > 1.5), there- 
fore the 0.267 arcmin^ of the Centre are expected to con- 
tain just ~ 7 red clump stars coming from the LMC field, 
which is far less than necessary to explain any of the features 
of its CMD. Indeed, the red clump in the Centre is made 
of 117 stars, which can be separated into the 89 "bright" 
ones {F814W < 18.75) - which correspond to the classi- 
cal red clump made of stars which likely passed through 
electron degeneracy in their cores - and 28 faint or "sec- 
ondary" ones (F814:W > 18.75) - which were likely able to 
avoid it. We can conclude, in a way similar to the case of 
NGC 419 (|Girardi et al.ll2009l . [ioTol ). that the probability 
that the dual red clump observed in the centre of NGC 1751 
is caused by LMC field stars is less than P ^ 5 x 10~^, and 
therefore negligible. 

Note that differences are quite evident in the position 
of CMD features between the Centre and Field, which are 
obviously mixed in the CMD of the Ring. The Field presents 
an old main sequence turn-off and subgiant sequence extend- 
ing to magnitudes as faint as FS14W < 21.5, and a younger 
main sequence extending as bright as F814VK < 16.5. Just 
traces of these features are present in the CMD of the Cen- 
tre. Moreover, the red clump in the Field do also present 
a composite structure, with a ratio of faint/bright stars of 
59/130. This latter feature is just expected from a field made 
of stars cover i ng a wide ran^e of ages and initial masses (see 
lGirardil[T999l : IPiatti et al.lll999l ). and being observed with 
very small photometric errors as in this case. 

2.2 Assessing photometric errors and 
completeness 

In order to characterize the errors in the photometry and 
the completeness of the sample, we have performed a series 
of artificial star tests (AST) on the reduced images (see e.g. 
iGallart et al .11 19991 : [Harris Zaritskvll200lh . 

The procedure consists of adding stars of known mag- 
nitude and colour at random places in each exposure, and 
redoing the photometry exactly in the same way as described 
in Sect. 12.11 The artificial stars are considered to be recov- 
ered if the input and output positions are closer than 0.5 
pixels, and flux differences are less than 0.5 mag. 

In order to avoid the introduction of additional crowd- 
ing in the images, artificial stars are positioned at distances 
much higher than their PSF width. So, our AST are dis- 
tributed on a grid spaced by 20 pix, which is each time 
randomly displaced over each set of exposures. 

A total of 1.04 X 10^ ASTs were performed, covering in 
an almost uniform way the CMD area of the observed stars 
as well as the area for which we build the "partial models" to 
be used in the SFH analysis (see Sect. 14.21 below). Then, the 
ratio between recovered and input stars gives origin to the 
completeness map of Fig. [H Note that the 90 % completeness 
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Figure 4. Completeness map, derived from the complete set 
of ASTs realised over NGC 1751 (centre plus ring areas), for 
both the F814W vs. F435W - F814W (left panel) and F814W 
vs. F555W-F814W (right panel). 
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Figure 5. Map of photometric errors as a function of input 
F435W, F555W and F814W (from top to bottom), as derived 
from the ASTs over the core of the cluster area (that is, in the 
Centre plus Ring). The errors are defined as the difference be- 
tween the recovered and input magnitudes. 



limit is located at F814W^ 24.5, which is well below the 
position of the MMSTOs in NGC 1751. 

Figure [5] illustrates the differences between the recov- 
ered and input magnitudes of the ASTs, as a function of 
their input magnitudes. These differences give a good han- 
dle of the photometric errors effectively present in the data. 
The error distributions are slightly asymmetric because of 
crowding. 




Figure 6. map for the Field best-fitting solution, as a function 
of distance modulus and V-band extinction. The continuous lines 
show the 68 % (black) and 95 % (white) confidence levels for 
the overall best-fitting solution, which is located at (m — M)o = 
18.50, Ay = 0.525. 



3 THE SFH OF THE LMC FIELD 
3.1 Overview of the method 

To recover the SFH from the ACS data, we use the pipeline 
built to analyse data fr om the VISTA su rvey of the Magel- 
lanic Clouds (VMC; see lCioni et"llll201lh. The method has 
been fully described and tested bv lKerber et al. (l2009l) using 
simulated near-infrared data, and by Rubele et al.l (|201Qr ) 
using ACS/HRC data for the SMC cluster NGC 419. The 
method consists in (1) building the Hess diagram for the 
data and a series of "partial models" which represent pop- 
ulations in limited interv als of age and metallicity, and (2) 
using the StarFISH code (^Harris Zaritskvl l200l[ l2004l ) to 
find the linear combination of partial raodels t hat minimizes 
a x^-like statistic as defined in lDolphinI (|2002l V The solution 
is characterized by the minimum x^, Xmin? and by a set of 
partial model coefficients corresponding to the several age 
bins. These latter translate directly into the star formation 
rate as function of age, SFR(t). 

The age-met allicity space occupied by the partial mod- 
els depends on the object under consideration. In the present 
work, we have to consider two distinct cases, corresponding 
to the cluster and LMC field as previously defined. This 
section deals with the field only. 

The LMC field is expected to follow a marked age- 
metallicity relation (AMR). This AMR has been measured 
by several auth ors using both field stars and star clusters 
of several a ges (Mackev G ilmordlioosl : iKerber et al.ll2007l : 
ICrocholski et al. . 2006 . 20071 In addition to the mean AMR, 
it is reasonable to expect a modest spread in metallicity at 
any given ag e. For this work, we adopt the scheme set by 
iKerber et al.l (|2009l l in which we build partial models at 



5 different metallicities disposed around the mean AMR of 
the LMC. Each partial model covers a range in logarithm of 
age of width 0.2 dex. For stellar populations younger than 
10^ yr, the numbers of observed stars are very small and 
hence we assume broader age bins, of widths 0.3 dex for 
7.2 ^ log(t/yr) ^ 8.0, together with a single age bin of 
width 0.8 dex for log(t/yr) < 7.2. The [Fe/H] separation 
between partial models is of 0.1 dex. 

For the initial mass function (IMF) we adopt the 
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Figure 8. Top panel: best-fitting SFH for the field, together with 
the random errors (Icr). Bottom panel: the mean age-metallicity 
relation. 



IChabrierl (|200lh one. The binary fraction is set to a value 
of 0.3 for binaries with mass ratios in the range between 0.7 
and 1.0, which is consistent with the prescriptions for bi- 
naries commonly used in works devo ted to recover the field 



SFH in the Magellanic Clouds (e.g. iHoltzman et al.[ll999l : 
iHarris Zaritskvll200ll : Ijaviel et al.il2005l : lNoel et al.ll2009h . 
Notice that this assumption is also in agreement with the 
few deterr ninations of binary fraction for ste llar clusters in 
the LMC (lElson et al.lll998al : IHu eralll2008l . both for NGC 
1818, a stellar cluster younger than ~ 100 Myr). 



3.2 The best-fitting solution 

Once the database of partial models is built, we run 
StarFISH to find the best-fitting solution to the observed 
CMDs, for a given value of distance modulus (m — M)o 
and extinction Ay. Both F814W x F435W - F814W and 
F814W X F555W — F814W Hess diagrams are used simulta- 
neously in the process of minimization. 

(m — M)o and Ay are then varied over the possible 
range of values. The map of Fig. [6] shows the results 
in the (m — M)o x Ay plane. The overall best-fitting solu- 
tion, with a Xmin of 1-6, is located at (m — M)o = 18.50 
and Ay — 0.525. The 68 % confidence level for this solution 
spans a narrow range in distance and reddening, which is 
just A(m — M)o = 0.12 mag and /\Ay = 0.07 mag wide. 

Figure [7] compares the Hess diagrams of the field data 
(left panel) and its overall best-fitting model (right panel). 
It is evident that the solutions found by StarFISH reproduce 
well the observed CMD features but for the Poisson noise. 

Finally, Fig. [8] presents the SFR(t) and age- 
metallicity relation (AMR) corresponding to this 
best-fitting solution. It is remarkable that the recov- 
ered SFR(t) presents features that are consistent with 
those commonly found in previous works, based on 



ISmecker-Hane et"aL 2002|: Javiel et al . 2005) and ground- 
based data ([Harris Zaritskvl l2Q0l[ 12009|). There is an 
initial burst of star formation followed by a quiescent period, 
with a marked and peaked star formation for ages younger 
than 4 Gyr (log(t/yr) = 9.6). In particular the peaks of 
star formation at approximately 1.5 Gyr (log(t/yr) = 9.2), 
500 Myr (log(t/yr) = 8.7), 100 Myr (log(t/yr) = 8.0) and 
10 Myr ( log(t/yr) = 7.0) are in ti ght agreement with those 
found bv lHarris Zaritskvl (l2009l ) for the global SFH of the 
LMC. Concerning the AMR, the result for the NGC 1751 
field is c onsistent with those derived from the LMC s tellar 
clusters (iKerber et all l2007l: [Harris Zaritskvl l2009l ) and 
for the LMC field (|Carrera et al.ll2008l ). 



4 THE SFH FOR NGC 1751 

4.1 Overview of NGC 1751 parameters from 
literature 

As for the LMC field, also for the NGC 1751 cluster we need 
a set of physical parameters to start with the SFH-recovery 
work. They are based on the following works: 

The cluster metallicity as de termined by the Ca II 
method is of [Fe/H] = -0.44 ±0.05 (ICrocholski et al.ll2006l l. 
which is a typical value for an intermediate- age LMC cluster 
( [Fe /H] = -0.48 d i 0.09; ICrocholski et alll2006l ) . 

iMilone et al] (|2009l ) identified a double MSTO in the 
HST/ACS F435W vs. F435W - F814W CMD for this clus- 
ter. Using isochrone fitting, they determined ages between 
1.3 and 1.5 Gyr, a distance modulus of 18.45 mag, Es-y = 
0.22 (Ay 0.70), and a metallicity of Z = 0.008 ([Fe/H] ^ 

-0.38)^ 

iMilone et al.l (|2009l ) also determine a binary fraction /b 
of 0.13 di 0.1 for mass ratios q larger than 0.6 in NGC 1751. 
Despite the small error bar quoted by them, their estimate 
is ad mitedly a prelimina ry one. The careful determination 
from lElson et al.1 (ll998b') for the LMC cluster NGC 1818, 
finds /b values varying from ^ 0.20 to ^ 0.35 as one goes 
from the cluster center to the outer parts. We adopt here 
the conservative value of /b = 0.2 for q > 0.7. Our previous 
results for NGC 419 (jRubele et al.ll2010l ) indicate that the 
results of the SFH recovery do not depend significantly on 
the choice of binary fraction. 

As for the extinction, the reddening maps from the Mag- 
ellan ic Clo uds Photometric S urvey (MCPS; IZaritskv et al.l 
20041) and IPeicha Stanekl (|2009l ) provide discrepant val- 
ues for the NGC 1751 direction. Within 3 arcmin from 
the cluster, MCPS gives Ay = 0.47 dz 0.34 for hot stars, 
and Ay — 0.59 di .39 f or cool stars. From the same 
dataset, IPessev et al.l (120081) determined Ay = 0.65 =b 0.06. 
Pejcha & Stanek (2009) instead find {Ey - i) = 0.150 ± 
0.293 {{Ay) = 0.293 ± 0.062), although their value is based 
on just 5 stars within a 2° x 2° area. 

The distance modulus to the LMC disk in the NGC 1751 
direction is of about 18 .55 mag, as revealed by indepen- 
dent methods: I van der M arel & Cioni ( 20Qll. AGB stars): 
Olse n Sal^ (|200l red clump stars); [Nikolaev et al.1 



f|200j . Cepheid stars') : Isiibramanian &: Subramaniam ( 201 c 
RC stars from MCPS). 

The above-mentioned works provide the initial guesses 
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Figure 7. The Hess diagram for the NGC 1751 Field as derived from the ACS data (left panels), as recovered by the best-fitting solution 
(central panels), and the map of residuals (right panels). 



for the many cluster parameters to be determined in this 
work. 



4.2 The partial models for NGC 1751 

For NGC 1751 we assume a constant age-metallicity re- 
lation, i.e., a single value for the metallicity for all ages, 
since so far there are no evidences for significant spreads in 
metallicity in such star clusters (e.g. Rlu cciarelli et al. 2008; 
iRubele et al.|[2010l V The age interval covered by our partial 
models goes from log(t/yr) = 8.9 to 9.4, which is much wider 
than the interval suggested by the position of NGC 1751 
MMSTOs. So, for each set of parameters, we have a total of 
10 partial models, completely encompassing the age inter- 
val of interest. We have explored 5 metallicity values, going 
from [Fe/H] from -0.75 to -0.35 at steps of 0.1 dex. For 
each one of these mean [Fe/H] values, a small metallicity 
spread of 0.02 dex is assumed. 



This definition of partial models might already be good 
enough to our aims to find the best-fitting solution for the 
cluster centre. However, we know that every portion of the 
ACS/WFC image is contaminated from the LMC field, and 
that this field contamination is well evident in the observed 
CMDs (especially for the Ring). Therefore, it is quite tempt- 
ing to add, to the above-mentioned list of partial models, an 
additional one corresponding to the LMC field. The hope 
is that this partial model will allow StarFISH to properly 
fit the field contamination across the CMDs on NGC 1751, 
hence improving the fitting of cluster itself. 

The inclusion of a partial model for the field is a nov- 
elty of this work, and is suggested as an alternative to the 
com monly used method of field star decontam ination (see 
e.g., iMilone et aLll2009l : iBonatto Bicall2007l ). which con- 
sists in subtracting from the cluster CMD the stars with 
colors and magnitudes similar to the ones in the field be- 
fore deriving the cluster parameters. The advantage of our 



The SFH of NGC 1751 7 



method is that the field becomes an integral part of the 
and error analysis; the latter is performed using the correct 
number statistics - i.e. considering the Poisson noise from 
the field, which is certainly present in cluster data - with- 
out implying any change in the method already set for these 
tasks. 

There are then two different ways at which this partial 
model for the field can be built. (1) The simplest one is that 
of taking the observed Hess diagrams for the field region 
(left panels in Fig. [Tj). This diagram however is affected by 
the Poisson fluctuations in the numbers of stars, so that 
it might not describe in a suitable way the field actually 
observed in other parts of the ACS images. (2) The second 
alternative is to use the Hess diagram built from the best- 
fitting solution of the field (right panels in Fig. [Tj) which 
is obviously much more continuous and smooth over the 
CMDs than the observed one. This model has another clear 
advantage: the Hess diagram can be easily re-built using the 
output SFH for the field together with the ASTs derived 
for the cluster Centre or Ring. In this way, we are able to 
simulate the field under the same conditions of crowding met 
in the cluster data. We indeed adopt this latter alternative 
in the following. 



4.3 The SFH for the cluster Centre 

4.3.1 Results with and without the LMC field 

The SFH-recovery is performed assuming a given set of 
(m — M)o, Av, and [Fe/H] values and fixing the binary frac- 
tion at a value of 0.2 in the case of cluster models. In order 
to limit the space of parameters to be covered, the proce- 
dure is essentially the following: for a given value of [Fe/H], 
we perform SFH-recovery for each point in a grid covering 
a significant region of the (m — M)o versus Ay plane, so as 
to build a map of the Xmin for the solutions. Examples of 
these maps are presented in Fig. [9l The maps are extended 
enough so that the minimum Xmin for ^ given value of [Fe/H] 
can be clearly identified, as well as the regions in which Xmin 
increases by a factor of about 1.5. The typical resolution of 
such maps is of 0.02 mag in (m — M)o and 0.02 mag in Ay. 

Let us first start discussing the case of the cluster Cen- 
tre. Fig.llOlshows the maps of Xmin cts a function of (m — M)o 
and Av, for two series of SFH-recovery experiments made 
under very similar conditions, i.e. using the same data and 
ASTs. The only difference is that in some cases (hereafter 
case A), we do not use the partial model for the LMC 
field in StarFISH, whereas in other cases we do it (hereafter 
case B). The result of considering the LMC field is quite 
evident: although in both A and B cases the best-fitting solu- 
tion is found for about the same value of (m — M)o and Ay, 
in case B the Xmin values are systematically smaller, which 
means better overall fits of the CMDs. Moreover, it is evi- 
dent that in case A the presence of the LMC field falses the 
determination of the best-fitting cluster metallicity: indeed, 
in case A the best-fitting model of [Fe/H] = —0.64 is found 
to be slightly favoured over the one with [Fe/H] = —0.44. In 
case B, instead, the best-fitting solution at [Fe/H] = —0.44 
is clearly favoured. Notice that, at the 1.5 Gyr old ages of 
NGC 1751, the field is found to have a metallicity of about 
—0.65 (see Fig. [8]), which probably helps, in case A, to move 
the Xmin minimum to [Fe/H] = —0.64. 



These experiments demonstrate that even a small frac- 
tion of field contamination may affect significantly the re- 
sults of SFH-recovery, if not properly taken into account. In 
the following, we adopt case B as the default, since it demon- 
strately takes the LMC field into account and improves the 
quality of the final results for the Centre of NGC 1751. 

4.3.2 Characteristics of the best- fitting solution 

Complete maps of Xmin for the Centre, as a function of 
(m — M)o, Av and metallicity, are presented in the left pan- 
els of Fig. [9l It may be noticed that the best solution is in- 
deed for [Fe/H] = -0.44, (m-M)o = 18.58, and Ay = 0.50, 
with a Xmin — 0.62. Such a small Xmin is already an indica- 
tion of an excellent fit to the observational data. 

This best-fitting solution and map of residuals are also 
presented in the Hess diagrams of Fig. 1111 Finally, the best- 
fitting solution for the cluster Centre is in the left panel of 

Fig.im 

4.3.3 Evaluating the errors 

To evaluate the errors for all involved parameters, the first 
step is to find the correspondence between the Xmin value 
for each model and its significance (or confidence) level, a. 
This correspondence was estimated simulating 100 synthetic 
CMDs generated with a number of stars equal to the ob- 
served CMD, using the best-fitting SFR(t) and its param- 
eters as the input for the simulations. So, after recovering 
the SFH for this sample of synthetic CMDs, it was possible 
to build the Xmin distribution and to establish the relation 
between the Xmin difference above the minimum and a. 

In the Xmin maps of Fig. [O] we superimposed the 68 % 
and 95 % significance levels for all the solutions for the 
centre. Only for the [Fe/H] = —0.44 map we find ample 
areas of the Ay versus (m — M)o diagram with solutions 
within the 68 % significance level of the best solution. Based 
on this figure, we determine (m — M)o = 18.58 d= 0.07 and 
Av = 0.50=1=0.05 for the cluster Centre (with random errors 
at the 68 % significance level). 

The left panel of Fig.[T2lshows the SFR(t) for the cluster 
Centre together with error bars. The most basic feature in 
this plot is that the SFR(t) is clearly non-null for three age 
bins, spanning the log(t/yr) interval from 9.05 to 9.2 (ages 
from 1.12 to 1.58 Gyr). Note that this result is not only valid 
for the best fitting model, but also across the entire 68 % 
significance level area of the Ay versus (m — M)o diagram. 
Moreover, it is non-null even in the case we adopt more 
restrictive limits for the random errors, i.e. if we plot the 
random errors for the 95 % significance level. 

Then, one may wonder how the solution for the Centre 
changes if we adopt a better age resolution in the SFH- 
recovery. This is shown in Fig. [TS] where we compare the 
solution for Alog(age) = 0.05 dex with the one obtained 
with the same data and methods, but for an age resolution of 
A log(age) = 0.025 dex. As we see, within the error bars the 
two solutions are essentially the same. The finer resolution 
in age is compensated by an increase in random errors. 

Therefore, we find evidence that in the NGC 1751 Cen- 
tre the SFR(t) has lasted for a timespan of 460 Myr. This 
is about twice longer than the ^ 200 Myr estimated by 
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Figure 9. Maps of the Xmin obtained from the SFH-recovery, as a function of (m — M)o and Ay, for several [Fe/H] values (from —0.34 
to —0.74 at steps of —0.1 dex, from top to bottom) and for both the cluster Centre and Ring (left and righ panels, respectively). The 
black lines delimit the regions within a 68 % (continuous line) and 95 % confidence levels (dotted lines) of the absolute best solution, 
which is found at —0.34 dex for the Centre, and at —0.74 dex for the Ring. The Xmin Centre best solution is of 0.625. 



iMilone et alj (|2009l l for the same cluster, based on a simple 
comparison with the MSTO locations of different isochrones. 

4.4 The SFH for the cluster Ring 

We have performed the same experiments of SFH-recovery 
separately for the Ring, as illustrated in the Xmin maps at 



the central and right panels of Fig. (9] In these cases, the 
levels of Xmin are significantly higher than for the cluster 
Centre. This result may seem surprising, considering that 
the Ring has a lower level of crowding and hence deeper 
photometry than the Centre. We consider as unlikely that 
these higher Xmin for the Ring could be simply caused by 
its higher level of contamination from the LMC field, since 
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Figure 10. Maps of the Xmin obtained during SFH-recovery in the Centre region, as a function of (m — M)o and Ay. The left panel 
shows the map for the best-fitting metaUicity of [Fe/H] = —0.64, obtained in the case A (i.e. not taking into account a partial model for 
the LMC field). The minimum X^i^i '^^ 0.77. The middle panel shows the same for case B (i.e. using the LMC field partial model) 
and [Fe/H] = —0.44, which is the best-fitting metallicity in this case. It is evident that these solutions are characterized by a significantly 
smaller level of X^\ri overall, with a minimum at 0.62. For comparison, the right panel shows the best-fitting solutions for case A and 
[Fe/H] = —0.44. Also in this case, the X^i-a significantly higher (and very close to the one in the letfmost panel). 




Figure 11. The Hess diagrams for the NGC 1751 Centre data (left panels), its best-fitting solution model (central panels), and the x^ 
map (right panels). The top panels are for the F435W - F814W vs. F814W diagrams, the bottom ones for F555W - F814W vs. F814W. 
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Figure 12. The blue lines show the SFR(f) for the cluster Centre (left panel) and for the Ring (right panel). The error bars are random 
errors. The green lines indicate the systematic errors, inside the area of 68 % confidence level in Fig. [O] 
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Figure 13. Comparison between the SFR(t) for two best-fitting 
solutions for the cluster Centre, obtained with different age reso- 
lutions: with Alogt = 0.050 dex (blue continuous line) and with 
Alogt = 0.025 dex (red dashed line). 



this field is very well modeled anyway, and fully taken into 
account in the SFH recovery. 

Instead, the main reason for the worst fits could be on 
the presence of a differential reddening of a bout Eb-v ^ 
0.10 m a^ within the ACS field, as found by iMilone et al.l 
(|2009l l and Goudfrooij et al. (in preparation). 

To test for the presence of differential reddening we have 



follow ed a similar procedure as described by iMilone et al.l 
following the position of fiducial lines in the 
F435W - F555W vs. F555W - F814W diagram, and as il- 
lustrated in Fig. 1141 So by means of the relative shifts in 
the fiducial lines in this colour-colour diagram along the 
reddening arrow we have found an extra reddening in the 
Ring region in relation to the Centre, with a magnitude in 
^F555 W - F814W similar to the one presented bv IMilone et al.l 
(l2QQ9l ). This extra reddening occurs prevalently in the bot- 
tom and upper extremities of the Ring in Fig. [1] 

Notice also that the extinction values for the Ring, 
found during the experiments of SFH-recovery are systemat- 
ically higher than the ones found for the Centre region (see 
Fig. O, which independently confirms the presence of dif - 
ferential reddening detected by us and IMilone et al.l (|2009h 
using fiducial lines. 

Unfortunately the recovered best-fitting (m — M)o and 
[Fe/H] values for the Ring are obviously spurious, since they 
are not consistent with the ones found for the Centre - which 
is a such small region that one can consider it free from 
differential reddening. Considering the high quality of the 
SFH-recovery for the cluster Centre (as demonstrated by the 
much smaller Xmin)? we assume that the correct distance and 
reddening of the Ring are the same ones as for the Centre, 
namely (m — M)o = 18.58 and Ay = 0.50. The right panel 
of Fig. [T2l shows the Ring SFH for these parameters. 

It is interesting to note that the Ring SFR(t), consider- 
ing just the random errors, seems to be slightly more spread 
in age than the one for the Centre. Indeed, the SFR(t) is 
found to be non-null in an additional, older age bin, as com- 
pared to the cluster Centre. However, when one considers 
the systematics errors in this bin, it is clear that this result 
is not solid. It becomes then impossible to take any conclu- 
sion in relation to a possible dependence of the SFR(t) with 
the cluster radius. 
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Figure 14. Colour-colour diagrams for the northern (left panel) and southern (right panel) half-sections of the Ring, illustrating the 
method used to determine the presence of differential reddening. Only the data for the lower main sequence, with F814W between 
21.5 and 23.5, is plotted (red dots). The blue dots mark the median position of the observed stars in small colour bins (dotted lines) 
perpendicular to the reddened main sequence. The continuous dark line is an unreddened model zero-age main sequence. The dashed 
line shows the position of the reddened main sequence that best fit these blue dots. Although the reddening vector (arrow) runs almost 
parallel to the main sequence, it is clear that the median position of the points in this plot can be used to derive the mean reddening in 
a field. We do it differentially, determining that between the northern and southern sections of the Ring there is a difference in Ay of 
0.30 mag, which corresponds to 0.10 mag in Eb-v- 



5 CONCLUDING REMARKS 

In this work, we confirm that NGC 1751 hosts a dual 
red clump, si milar to the one ob served in the SMC clus- 
ter NGC 419 (jGirardi et alJ 120091 ). Then, we perform SFH- 
recovery via the classical method of CMD reconstruction 
with the sum of single-burst stellar populations. The basic 
result is that in the NGC 1751 Centre star formation is found 
to last for a timespan of 460 Myr. The same result is con- 
sistently found irrespective of the method used to take field 
contamination into account, of the age resolution adopted, 
and for a significant region of the Ay versus (m — M)o plane. 
The age resolution of the method is at least 3 times smaller 
than this interval. In addition, the best-fitting model is in- 
deed a quite good representation of the observed data. On 
the other hand, the CMD for the cluster Centre is so obvi- 
ously sharp and clean from the LMC field, that it is hard to 
imagine that some important effect has not been properly 
considered in our analysis. 

Our re s ults contrast with those obtained by 
iMilone et al.l (|2009h , who identify just two episodes of 
star formation separated by 200 Myr in this cluster, using 
the method of isochrone fitting. 

For the cluster Ring, the results indicate a SFR times- 
pan of at least 460 Myr, with a significantly lower quality 
of the results, as indicated by the larger Xmin sind by the 
best fitting model being found for a distance different from 
the cluster Centre. These failures probably depend on the 
presence of significant differential reddening across the Ring. 
So, we prefer not to take any conclusion from this region. 
Unfortunately, our method finishes in not providing any in- 
dication about possible variations in the spatial distribution 
of the stars with different ages, which would be important 



for the goals of understand ing the m echanism of extended 
star formation fsee Goudfr ooii et al.|[^009,i V 

Together with our previous findings for NGC 419 
(jRubele et al.1 l2010l ). the results for the NGC 1751 Cen- 
tre argue in favour of multiple star formation episodes 
(or continued star formation) being at the origin of the 
MMSTOs in Magellanic Cloud clusters with ages around 
1.5 Gyr. The hypothesis of a spread caused by the pres- 
ence of fast rotators among the upper main sequence 
stars ([Bastian de Mink 2009 ) is se parately discussed by 
iGirardi, Eggenberger, k, Miglio (|201lh , who conclude that it 
does not offer a valid alternative to the conclusions reached 
in the present paper. 
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